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Abstract 

We study the Cu-0 valence instability (VI) and the related phase separa- 
tion (PS) driven by Cu-0 nearest-neighbor repulsion U p d, using an effective 
extended one-band Hubbard model (H e ff) obtained from the extended three- 
band Hubbard model, through an appropriate low-energy reduction. H e ff is 
solved by exact diagonalization of a square cluster with 10 unit cells and also 
within a slave-boson mean-field theory. Its parameters depend on doping for 
Upd / or on-site O repulsion U p ^ 0. The results using both techniques 
coincide in that there is neither VI nor PS for doping levels x < 0.5 if U p d < 2 
eV. The PS region begins for U p d > 2 eV at large doping x > 0.6 and in- 
creases with increasing U p d- The PS also increases with increasing on-site Cu 
repulsion Ud- 

PACS numbers: 74.72.-h, 74.25.Jb 
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I. INTRODUCTION 



It is widely accepted that the electronic properties of the cuprates are well described by 
the three-band Hubbard model H 3b |l]-|3| . Zhang and Rice 0] suggested that this model can 
be reduced to an effective one-band t — J model under certain simplifying hypothesis: small 
O-Cu hopping t pd , and zero on-site O Coulomb repulsion U p , Cu-0 Coulomb nearest-neighbor 
interatomic repulsion U p d and 0-0 hopping t pp . While the one-band effective models can 
explain the magnetic properties of high-T c materials and can provide magnetic mechanisms 
for superconductivity these models, with fixed parameters, do not take into account 

the effects of the Cu-0 repulsion U p d- The density of carriers in cuprate superconductors 
is very low. The average distance between two holes is larger than 7 A and the Cu-0 
distance is ~ 1.9 A and therefore U p d is expected to be poorly screened. Neglecting screening 
U p d ~ 7 eV, while direct calculations give U p d ~ 3 eV |7j and constrained-density-functional 
results predict U p d ~ 1 eV ||. Thus, as pointed out early after the discovery of high-T c 
superconductivity, the effects of U p d can be important @||. These effects include a charge- 
transfer mechanism for superconductivity [3fl|10| , and marginal -Fermi-liquid behavior |TT 



It is known that a sufficiently large U p d induces a charge-transfer instability (CTI) and 



a valence instability (VI) related with Cu-0 charge transfer |I2" |. These instabilities are 
associated with phase separation because of the coupling between the fluctuations of the 



valence and the total density fllqJLl 



The three-band Hubbard model including U p d has been studied by exact diagonalization 
of small clusters ||, random-phase approximation |3| JT3| , p!5|Jl6| , Gutzwiller variational wave 



function [17], and slave-bosons with 1/N expansion, for the particular case Ud = oo and 



with U p and U p d treated in the Hartree approximation |14],[T8],[19| . The results of these works 
are qualitatively similar in general. Increasing U p d the system reaches the CTI, VI and the 
related phase separation. Near the phase-separation boundary, the effective interaction in 
the Cooper s and d channels becomes attractive. 

Bang et al. in weak coupling |13|] and Raimondi et al. in strong coupling JE1J have 
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shown that even though the energy uo exc of the charge-transfer collective excitonic mode 
at zero wave vector, decreases with increasing U p d, it remains finite at the CTI. However, 
this mode is coupled to the zero-sound mode, leading to a charge-transfer-mode mediated 
attraction which allows the violation of the Landau stability criterion Fq > — 1. This leads 
to the simultaneous divergence of the compressibility and charge-transfer static susceptibility 
indicating that phase separation and the CTI take place. 

In spite of the qualitative agreement on the facts mentioned above among the differ- 
ent techniques, the regions of phase separation differ in the strong and weak coupling ap- 
proaches, and the critical value of U p d where the instabilities take place depends strongly 
on the approximation used. Also, the effect of the variation of the parameters has not been 
investigated in detail and realistic values of Ud have not been studied so far. Thus, further 
research on these subjects seems necessary. 

On the other hand, it is important to address the question of to what extent the above 
mentioned properties of the three-band Hubbard model can be accounted for using a one- 
band effective model derived from the former through a sufficiently accurate mapping proce- 
dure. Generalized t — J and one-band Hubbard models have been obtained performing sys- 
tematic low-energy reductions of the three-band Hubbard [E0-E7). These reductions either 



are based on a cell method [p0| , plj| or use an effective spin-fermion model with renormalized 



parameters as an intermediate step in the derivation p7| , [28[| . In contrast to the original 
derivation of Zhang and Rice, these methods allow to extend the mapping procedure to 
realistic and large values of Cu-0 hopping t p d, and show that the Zhang-Rice singlets are 
stable for large Cu-0 covalency. There have been an important amount of research devoted 
to the study of the validity of the low-energy reduction and the stability of the Zhang-Rice 
singlets. In particular, the roles of local triplets states [^9| and apical O ions []30| have 
been recently investigated. Exact diagonalization of small clusters have shown that the low- 



energy spectrum 0,[n|] and magnetic properties |32|] of the three-band model can be well 
reproduced by a t — t' — t" — J model, where t' is the next-nearest-neighbor hopping and t" 
is a next-nearest-neighbor hopping combined with a nearest-neighbor spin flip. This three- 



body term, with the sign corresponding to large 0-0 hopping t pp in the original three-band 
model, can stabilize a resonance- valence-bond superconducting ground state |J. 

In this paper, starting from the three-band Hubbard model, we derive an effective gen- 
eralized one-band Hubbard model , using the cell method generalized to take into account 
properly the intercell part of the O intratomic repulsion U p and the Cu-0 interatomic re- 
pulsion Upd- As a consequence of these interactions, the parameters of the effective model 
become dependent on the particle density. Solving the effective one-band model exactly in a 
VTO x vTO cluster, and also in the slave-boson mean-field theory of Kotliar and Ruckenstein 
f33fl we study the above mentioned valence instability (VI), the regions of phase separation 
(PS), and also the dependence of the effective one-band parameters with U p d and doping. 
The advantage of our approach in comparison with other methods discussed above for the 
study of the VI, CTI and PS, is that the intracell correlations, including Cu on-site repulsion 
Ud and a large part of the U p d and U p terms, are taken into account exactly at each cell. This 
allows us in particular to take realistic values of Ud (~ 10 eV), while in previous treatments 
Ud was either small or infinite. 

In section II we describe the different Hamiltonians and the method for low-energy re- 
duction. Section III contains the results and section IV is a short summary and discussion. 

II. REDUCTION FROM THE THREE-BAND TO A ONE-BAND HAMILTONIAN 

We start from the three-band extended Hubbard model: 

H 3b = A PjaPj<? + u d d % d iAi d H + Up^PjiPttPjiPjl 

j i j 

+ U pd d\a> d ia'P\ + 5^+5<r + tpdj^ipl+gji, + k.C.) ~ t pp ^ p] +ya Pja (1) 

iSao' iScr jja 

The sum over i (j) runs over all Cu (O) ions. The vector 5 (7) connects a Cu (O) site with 
one of its four nearest O atoms. The operator d\ a (Pj a ) creates a hole with symmetry d x 2_ y 2 
(pa) at site i (j) with spin a. The phases of half of the orbitals have been changed in such 
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a way that for all directions, the hoppings are positive (t p d, t pp > 0). The parameters of the 
model are known approximately from constrained-density-functional calculations [Q. 

The first step in the cell method [^1],[24],[26|] is to change the basis of the O orbitals to 
linear combinations which hybridize (oikv ) and do not hybridize (7^0- ) with dk a orbitals, 
due to the term in t pd at each point k of the reciprocal space. We denote the Wannier 
functions of these orbitals as and 7j CT respectively. The 7^ (non-bonding states) lie very 
high in energy and are neglected. After the change of basis, the original Hamiltonian can be 
separated in two parts, one containing the intracell (and generally larger) terms and another 
containing the intercell terms: 

-^36 = ^2 Hi + H inter , (2) 
% 

with 

Hi = (A- n(0)t pp ) ®L a iv 

(7 

+U d d\ ] d, iA d\ l d ii + 2t pd X(0) y^d^a^ + h.c.) 

(7 

+U pd f(0) 22 4a' d ^' a L a ^ + Uph^a^a^al^ii (3) 

era' 

Hinter = ^2[2t pd X(j - i)d j ja a ia - t pp fi(j - i)a] a a i<T + h.c] 

+ U pd Y [/(* -J^- l ) n t' a ]a a l<r + h.C.) 
I 

+U P ^2 h(i - j,i-l,i- m)a^a m] a} jL ai l (4) 

ijlm 

The functions of the lattice vectors A, //, v, f and h are given in Ref. They decay 

rapidly with increasing argument and as a consequence, most of the original hoppings and 
interactions are contained in ^ Hi. The ^ indicates that in the last sum the term with 
i = j = m = I is excluded. 

The ordinary cell method consists in solving Hi exactly in the subspaces of 0, 1 and 2 par- 
ticles and retaining the ground state in each subspace. The non trivial retained eigenstates 
have the form: 



\ta) = (B 1 dl-B 2 al)\0). (6) 

The energies of these states will be denoted E 2 and E\ respectively. The first term of |z2) 
corresponds to the Zhang- Rice singlet in the original derivation M. Our modification consists 
in considering the coefficients as variational parameters which are determined minimizing 
the total energy. We treat the terms in U p d of H in t e r with i / j / I replacing nf a , by 
its expectation value. Using the fact that ^ f(i — j,i — I) = 0, this is equivalent to the 
treatment of Ref. ||22| . Similarly in the terms in U p of Hi nter with i = m and j ^ I we 



replace o^ ia OL ia by its expectation value. The same treatment is done for i / m and j = I. 
The remaining terms except those with i = m and j = I were neglected. The intercell 
Hamiltonian is written as a function of the retained eigenstates of ifj, and mapping the 
latter into those of the one-band Hubbard model at site i, one obtains the following effective 
model: 

H = Ei 2j n ia + U 2j n^n^ + ^ {c) a c i<T {t AA (l - m - a ){l - nj 

ia i <ij>a 

+tAB[ni-a(l - + "i+J ,-<x(l ~ «i-<r)] + t B Bni- a rij _ CT + /l.C. 

^ Vii(l - rii-crXl - nj- a ')n i(T ri ja i 

<ij>aa' 

^ Vi 2 ((l - ni-ajrij-a* + rii- a (l - n jt _ a/ ))n ia n ja > 

<ij>oo' 

^ l^7^_ (T n 3 -_ (T /n io .njv (7) 

<ij>CT(T' 

We obtain the ground state energy of if as a function of the three free variational 
parameters (Eqs. || and || ) using two different approaches: 



1. The slave-boson approximation (SB) of Kotliar and Ruckenstein [p5| , |33f after treating 
the nearest-neighbor repulsions in the Hartree-Fock approximation. 

2. Exact diagonalization (ED) of a square cluster with 10 unit cells. 



III. RESULTS 



We begin this section by studying the parameters of the one-band effective model and its 
dependence with the original parameters and doping. As a basis for our study we take the 
parameters of the original three-band Hubbard model determined from constrained-density- 
functional approximation M. We take the unit of energy as t p d = 1 — 1.3 eV. However 
we also study the effect of U p d and Ud within a wider range. In the (unrealistic) limit U p 
= U p d = 0, the intercell interactions vanish (Vn = Vu = V22 = 0) and the variational 
parameters A4, Bi (Eqs. |5] and |5] ) which minimize the total energy correspond to the 
eigenstates of the cell Hamiltonian Hi (Eq. |]) for all dopings. Thus the effective one-band 
parameters are independent of doping. Increasing U p d and U p , the intercell interactions 
appear and the variational parameters start to differ from those corresponding to the low- 
energy eigenstates of Hi to take into account better the intercell interactions. The variational 
parameters as well as the one-band effective parameters become doping dependent, although 
this dependence is very weak for small U p d and U p . 

In Figs. 1-4 we show the one-band parameters as functions of doping for different values 
of U p d using ED. The effective on-site repulsion U increases with U p d except for small values 
of A and high doping where the amount of O states in the local singlet (Eq. ^) increases and 
the mixed Cu-0 part (which pays ~ U p d) decreases. Related with this fact, U decreases with 
doping for small values of A. This reflects a metallization of the system with a larger amount 
of Cu-0 covalency which should be reflected in an increase of the conductivity. This might 
be related to the sudden increase in the Hall conductivity observed in L^-zSr^CuC^ for 
x > 0.17 [Q. As a consequence of the larger amount of O holes in the doubly occupied cells, 
and the change in the variational parameters Bi (Eq. ^) to minimize the loss of energy due 
to the Cu-0 repulsion U p d, the energy E\ and the amount of O holes in the singly-occupied 
cells increase with doping for small A and large Upd- For others values of A and U p d, E\ is 
rather insensitive to doping (see Fig. 2). 

For small values of U p d the three correlated hoppings of the effective one-band model are 
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rather insensitive to doping and similar in magnitude. For large values of U p d the last two 
terms in Eq. (|]) can not be neglected and affect somewhat the magnitude of these hoppings. 
As for the case of the correlated hoppings, if U p d < 1 the different effective nearest-neighbor 
repulsions are similar between them and rather independent of hopping. However in contrast 
to the former, the V l s increase significantly with U p d- For small A and large U p d, in spite of 
the fact that there is a significant charge transfer from Cu to O with doping, the effective 
repulsions do not change very much because the decrease (increase) in Cu-0 repulsion is 
approximately compensated by the increase (decrease) in 0-0 repulsion. 

In Figs. 5 and 6 we show the dependence of the expectation value rid = (n^ + n ii) 
on A. As expected, rid increases with A. For small U p d the results obtained by exact 
diagonalization (ED) and slave bosons (SB) are qualitatively similar, although there are 
quantitative differences of the order of 20% for larger dopings. For larger values of U p d 
there are some important qualitative differences between both methods which are discussed 
below, however both methods coincide in that for intermediate and small values of A, rid 
decreases with doping as a consequence of Cu-0 charge transfer which overcompensates the 
effect of doping (see Figs. 6 and 7). This is related with the increase in Cu-0 covalency and 
the metallization mentioned before. In contrast, for large values of A, rid increases faster 
with doping than for small values of U p d- As a consequence of these opposite effects of U p d 
for small and large A, the dependence of rid with A becomes more abrupt for large U p d, 
eventually driving the valence instability (VI). The effect of reducing Ud induces a larger 
Cu occupancy at large values of A favoring the VI. The static valence susceptibility is 
defined by: 

_ d(n d - rip) | _ 9 dnd l 
Xv- dA \n-^ dA \n, \p) 

where n is the total occupation per cell and n p = n — nd is the O occupation per cell, xv 
differs from the charge-transfer susceptibility xct because the latter is defined at constant 
chemical potential instead of at constant n. The valence susceptibility xv is shown in Fig. 
8. For each value of x — n — 1, there is a critical value of U p d (US d (x)) for which xv diverges 
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at a critical value of A (A c (x)), indicating the presence of a VI. For U p d > U^ d there is a 
discontinuous transition as shown for example in Fig. 6 for x = 0.8. The minimum possible 
value of Up d (x) occurs for maximum doping (x — 1). The critical value U° d (x) increases with 
decreasing x. Note that U pd (x) found with SB is always lower than that found with ED. 
Also, at x = the SB approximation gives an artificially large \v & t the metal-insulator 
transition due to the vanishing of the double occupancy (d) in the insulating phase and 
the exaggerated large increase of d with A in the metallic phase near the metal-insulator 
boundary. Within the ED, d grows smoothly with A in absence of the VI as one expects. 

In Fig. 9 we show the phase-separation diagram calculated with the SB approximation. 
Contrary to previous strong-coupling mean-field approximations, where always one of the 



phases between which phase separation takes place has x = |14[ , we obtain, in qualitative 
agreement with weak-coupling approximations |TjJ that phase separation takes place be- 
tween a phase with large doping x = 1 and another one with x depending of the parameters. 
As expected, the region of phase separation grows with increasing U p d- At least for U p d < 1,( 
i.e. for parameters near those obtained by constrained-density-functional approximation ||) 
there is no phase separation. The effect of decreasing XJ d on phase separation is to suppress 
it for large values of A. This is related with the change of the character of the local singlet 
|i2) (Eq. HD from mainly Cu-0 or 0-0 to a doubly occupied Cu orbital. For low values of A, 
the phase diagram is rather independent of Ud- Also, the effects of U p d and Ud on the phase 
diagram on a qualitative level can be understood on general physical grounds: increasing 
these repulsions favors localization, inhibits the kinetic energy terms, and as a result the 
dependence of the energy on occupation is more flat, favoring phase separation. Although 
using ED we have only a few possible different densities, due to the small size of the cluster, 
the results for the compressibility at these densities agree qualitatively with the SB phase 
separation diagram. However there are some important differences in a quantitative level: 
a) there is no phase separation for U p d < 2. b) For U p d = 3 there is a small island inside 
which the system phase separates in two phases, both with compositions inside the interval 
< x < 1. Some of these features can be seen in Fig. 10, where we show the chemical 
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potential jj, as a function of doping. In table I we list the densities for which the compress- 
ibility is negative (a sufficient but not necessary condition for phase separation) according 
to our ED results for different parameters. 



IV. SUMMARY AND DISCUSSION 

In this work, we have studied the effects of Cu-0 Coulomb repulsion U p d on valence 
instabilities and phase separation. Using a variational form of the cell method, we derive 
a one-band effective model with doping dependent parameters. This model has the form 
of a generalized Hubbard with correlated hoppings and nearest-neighbor repulsions. This 
approach allows us to study realistic values of the on-site Cu repulsion Ud ~ 10 eV and is 
accurate enough for moderate values of U p d and on-site O repulsion U p . We obtain that 
for parameters near those derived using constrained-density- functional approximation ||, in 
particular U p d ;$ 1-3 eV, the low- energy physics of the three-band Hubbard model can be 
well represented by the effective one-band model, with doping independent parameters. The 
effect of U p d is merely to renormalize the effective parameters and there are neither valence 
instabilities nor phase separation. 

For larger U p d the effective parameters become strongly doping dependent and for U p d ~ 3 
eV valence instabilities and phase separation appear. In agreement with previous weak- 



coupling results flql , we obtain that phase separation starts to occur for large doping values 
(x ~ 0.7 ±0.1 according to the exact-diagonalization results or x ~ 1 according to the slave- 
boson ones). The valence instability begins at x — 1 in both treatments. This facts are 



in contrast to results of alternative previous strong-coupling approaches JT2|JTJ,[T8| . This is 
probably due to the facts that on the one hand, the slave-boson approximation produces an 
artificial increase of xv at the metal- insulator transition (x — 0), as discussed in section III, 
favoring PS at low x. On the other hand in Refs. p^fl^JT9|| , the U p d and U p terms were taken 
in a mean-field approximation, essentially equivalent to Hartree-Fock , (and since Ud is large, 
this induces an artificially large increase of the energy with doping), while in our present 
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treatment most of the U p d term and a large part of the U p one are treated exactly. In fact, 
we have verified that treating both terms in Hartree-Fock before performing the low-energy 
reduction to the one-band model (as was done in previous studies of the metal-insulator 
transition |[25||), the results become qualitatively similar to those of Refs. fT2|,|T4| . 
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FIGURE CAPTIONS 



Fig. 1: On-site repulsion of the effective one-band model (Eq. |7]) obtained by exact 
diagonalization (ED) (see text) as a function of hole doping x for different values of the 
Cu-0 charge transfer energy (A) and Cu-0 repulsion: U p d=l (solid), 2 (dashed), 3 (dotted), 
4 (dotted-dashed). Other parameters of the original three-band model (Eq. |^ ) are: U d = 7, 
Up 4, tpd 1, tpp .5. 

Fig. 2: One particle on-site energy of the effective one-band model (Eq. |^) as a function of 
doping for different values of A and U p d- The three lower curves correspond to A = 1 and 
the other three to A = 3. Other parameters as in Fig. 1. 

Fig. 3: Nearest-neighbor hoppings of the effective one-band model (Eq. |7|) as a function 
of doping for different values of A ((a) A = 1, (b) A = 3) and U p d- Other parameters as in 
Fig. 1. 

Fig. 4: Nearest-neighbor Coulomb repulsions of the effective one-band model (Eq. |7|) as 
a function of doping for different values of A and U p d=l (solid), and 3 (dotted). Other 
parameters as in Fig. 1. 

Fig. 5: Cu occupancy as a function of Cu-0 charge transfer energy (A) for different dopings 
(indicated at the right of the corresponding curve) obtained by (a) exact diagonalization and 
(b) slave bosons. Other parameters are U p d = 2, Ud = 7, U p = 4, t p d = 1, t pp = 0.5. 
Fig. 6: Same as Fig. 5 for U pd = 4, U d = 10. 

Fig. 7: Cu occupancy as a function of doping for different values of A (indicated at the 
right of the corresponding curve) and U d , obtained by slave bosons (SB). Other parameters 
are U pd = 3, U p = 4, t pd = 1, t pp = 0.5. 

Fig. 8: Valence susceptibility as a function of A for the different positive dopings of Fig. 6 
obtained by (a) exact diagonalization and (b) slave bosons. Other parameters as in Fig. 6. 
Fig. 9: Phase-separation diagram A vs. x obtained using SB for different values of U d and 
Upd- The dotted line indicates the first-order valence transition for U p d = 4 (for U pd = 2,3 
there are no first order VI transition). The solid squares represent the same obtained with 
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ED. Other parameters as in Figs. 1,5. 

Fig. 10: Chemical potential as a function of doping for several values of A calculated with 
(a) SB and (b) ED. Other parameters are U p d — 3, Ud — 7, U p = 4, t p d = 1, t pp = 0.5. 



VI. TABLE CAPTION 

Table I: Values of the density n for which the compressibility is negative according to ED 
for different values of Ud, U v d and A. Other parameters are U p = 4, t p d = 1, t pp = 0.5. 
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